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The theory and design of superbackscattering nanoparticle dimers are presented. We analytically 
derive the optimal configurations and the upper bound of their backscattering cross-sections. In 
particular, it is demonstrated that electrically small nanoparticle dimers can enhance the backscat¬ 
tering by a factor of 6.25 with respect to single dipolar particles. We demonstrate that optimal 
designs approaching this theoretical limit can be found by using a simple circuit model. The study 
of practical implementations based on plasmonic and high-permittivity particles reveal that fourfold 
enhancement factors might be attainable even with realistic losses. 


I. INTRODUCTION 

The scattering of electromagnetic fields by subwave¬ 
length particles is of fundamental interest for a wide 
range of disciplines including physics, optics and engi¬ 
neering [1-5]. In most cases, the response of electri¬ 
cally small particles can be approximated by the exci¬ 
tation of an electric dipole, usually aligned with the in¬ 
cident electric field, and whose value is directly propor¬ 
tional to the local electric field and particle volume [I, 2]. 
Despite this fact, researchers have aggressively pursued 
the excitation of more exotic and sophisticated scatter¬ 
ing responses. These include, among others, magnetic 
dipoles [6-11], higher order modes (HOMs) [12-16], as 
well as their interdependence via magnetoelectric cou¬ 
pling [17, 18]. When properly combined, the excitation 
of these more advanced responses enables unusually large 
scattering behaviors, including the so-called superscat- 
tering [19, 20], and superdirective scattering [21-23]. Ad¬ 
ditional degrees of freedom can be introduced by using 
nonlinear elements [24]. On the other hand, the total 
[25-27] or angularly selective [28, 29] suppression of scat¬ 
tering is equally important, as well as the scattering min¬ 
imization while maintaining a useful amount of absorbed 
power [30-32]. In addition to the research on specific 
configurations, there have also been substantial efforts 
in understanding the fundamental limits and physical 
bounds of both broadband [33-36] and time-harmonic 
[37-40] scattering by passive particles. Even if the limits 
of passive particles are too stringent, the scattering can 
be boosted by using active media, resulting in the devel¬ 
opment of subwavelength nanoparticle lasers [23, 41, 42]. 

Despite all of this active research in the field of scat¬ 
tering physics and engineering, little to no attention has 
been paid to enhancing the backscattering of subwave¬ 
length particles (i.e., to the reflection from, or to the 
power re-radiated against the direction of propagation of 
the incident wave). This fact is surprising in view of the 
broad scope of particles with large backscattering cross- 
sections. In essence, the operation of any spectroscopy, 
communication, remote sensing (radar, sonar, ...), ma¬ 
nipulation and/or imaging system that makes use of a 
single emitter/receiver device relies on a backscattering 


measurement. This wide range of applicability should 
undoubtedly motivate the research on particles featuring 
enhanced backscattering cross-sections. 

In a recent work, the authors have derived an upper 
bound for the time-harmonic backscattering cross-section 
(normalized to the wavelength squared) of a passive par¬ 
ticle [37] 

^ < — Dscat ^ -Dgcat (^) 

This bound illustrates the limit on backscattering as a 
function of the scattering directivity in the forward and 
backward directions, Dscat(ki) and Dscat(~k/, respec¬ 
tively (i.e., stands for the direction of the incident 
wave). Intuitively, in order to maximize the scattering 
along a given direction, in this case the backward direc¬ 
tion, the particle must not only focus its scattering along 
this direction, but also along the incident direction. This 
backward-forward requirement is necessary in order to 
produce the destructive interference required to extract 
energy from the incident field. Eollowing this philosophy, 
one can fully align the design of superbackscattering par¬ 
ticles with the superdirective beamforming methods de¬ 
veloped in antenna engineering. In fact, this concept has 
already been successfully applied to the design of super¬ 
backscattering antenna arrays at microwave frequencies 
[43]. Eollowing a similar philosophy, the present work in¬ 
vestigates the design of superbackscattering nanoparticle 
dimers. These geometries have been selected in view of 
previous positive experiences with two-element superdi¬ 
rective antenna arrays [44-47] and plasmonic dimers [48]. 

II. THEORY 

Let us then begin by examining a nanoparticle dimer 
composed of two nanoparticles symmetrically located 
with respect to the origin on the z-axis with inter-particle 
separation d (c.f., figure 1(a)). The dimer is illumi¬ 
nated by a time-harmonic time dependence as¬ 

sumed throughout) x-linearly polarized plane-wave prop¬ 
agating along the z-axis, i.e., E* = xEq exp(—j ko z). Eor 
subwavelength particles whose separation is larger than 
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FIG. 1. (a) Sketch of the scattering geometry: a plane wave 
illuminates a nanoparticle dimer, exciting the electric dipole 
moments of each nanoparticle, (b) Equivalent circuit model. 


the particle radius, i.e., d — (ri + r 2 ) > ri,r 2 , the fields 
scattered by the nanoparticles can be fairly approximated 
by those radiated by an electric dipole p, or, equiva¬ 
lently, an electric Hertzian dipole (EHD) with current 
moment Id\ = p/jcj. Under this dipolar approximation, 
the particle dimer can be treated as a two-element lin¬ 
ear array of EHDs. In the far zone, the field scattered 
by the dimer consists of a spherical wave with electric 
field W = Eof (r) exp(—j ko r )/ {k^ r), where f (?) is the 
(unit-less) scattering pattern of the dimer, given by 


( .knd^^ .knd^^\ 

( 2 ) 

where fdp is the dipolar scattering pattern of each indi¬ 
vidual particle, and Ai and A 2 are the excitations coef¬ 
ficients describing the relative magnitude and phase of 
the dipole excited at each nanoparticle. Consequently, 
the scattering directivity of the nanoparticle dimer can 
be written as follows 


-^scat (^^) 
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In view of (1), the least upper bound of the backscat- 
tering cross-section of a nanoparticle dimer can be 
found by identifying the optimal excitation coefficients, 
Ai = A^^^ and A 2 = ^ 2 ^^, that maximize the forward- 
backward directivity product, Unseat(~ki) ^scat(ki)- In 
general, different excitation coefficients are obtained for 
each separation distance d, i.e., A^^^ = A^^{d)^ for 
p = 1,2. In particular, it is found in Appendix A that, 
in the limit of electrically small separation distances 


{k^d 0), the optimal excitation coefficients to maxi¬ 
mize the backscattering are given by 

Ar{d)=-Ar{d) = ^^ ( 4 ) 


where Uq G C is an arbitrary complex constant. 

Interestingly, the superbackscattering response is ob¬ 
tained exactly at A^^^ = — i.e., the electric dipole 

moments excited in each nanoparticle must be of the 
same magnitude but opposite direction. This simple con¬ 
dition has been previously mentioned in the literature 
under other circumstances, and it deserves a careful dis¬ 
cussion. To begin, superdirective antenna articles relate 
opposing currents with radiation superdirectivity [44], 
which, in scattering terms, is directly associated with 
superdirective scattering in the forward direction. How¬ 
ever, the maximum of the forward scattering directivity, 
also reported in the Appendix A, is obtained for slightly 
different coefficients = —ao/ {k{)d) and = 

ao [1/ {kod) + j2/5]. Note that the coefficients for opti¬ 
mal backscattering and optimal forward-scattering be¬ 
come very similar for small separation distances, i.e., 
A{^d ^ A 2 ^d ^ A^^^ despite the fact that they 

lead to significantly different patterns and scattering 
cross-sections. Moreover, the fact that significantly dif¬ 
ferent responses are obtained with numerically similar co¬ 
efficients emphasizes the stringent fabrication tolerances 
associated with the construction of superdirective devices 

[45]. 

It is also worth noticing that, when condition (4) 
is satisfied, the net electric dipole moment of the sys¬ 
tem: Pi + P 2 , is zero. This effect is usually associated 
with scattering cancellation configurations in many body 
problems [49]. Intriguingly, our analysis remarks that the 
backscattering is optimized when the predominant scat¬ 
tering response of the dimer is canceled out. Naturally, 
this is only possible if the usually secondary responses are 
strengthened. In particular, note that opposed electric 
currents have been traditionally associated with strong 
magnetic dipole moments [50]. Other works have also 
emphasized the role of the electric quadrupole excitation 
in these structures [51, 52]. To shed more light onto how 
the secondary responses are responsible for the enhanced 
backscattering, it is worth examining the multipolar de¬ 
composition of the above structure. 

Without any loss of generality, the scattered field out¬ 
side the particle can be written in multipolar form as 
follows: 

E* (r) = ^ (r) + (r)] (5) 

{A 

A complete description of the above notation can be 
found in [40]. For the configuration depicted in fig¬ 
ure (1) (a), the only nonzero scattered field coefficients 
correspond to the set of even TM modes, and 

odd TE^ TE^^^^ modes. These coefficients can be writ- 
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ten as follows 

a^™ = [hdl - (- 1 )” hdl] {2 (n + 1 ) {kodf-^ ipn} 

(6) 

= [(-!)" hdl + hdl] { {kodf' tpn } (7) 

with the term 

_ 1 1 

^ 47r 2^^ n (n + 1) (2n — 1)!! 


III. CIRCUIT MODEL 

After identifying the magnitude and phase of the dipole 
moments that lead to a superbackscattering response, 
one must still find the correct design (nanoparticle geom¬ 
etry and materials) that provides the optimal excitation 
coefficients (4). We now demonstrate that this can be ac¬ 
complished by using a simple circuit model. Specifically, 
the local field equations acting on each particle can be 
written in circuital terms as: 

Zih = Eodle^’^ ^Zch ( 11 ) 


It is apparent from ( 6 ), (7) that the first dominant term 
is the electric dipole, which is oc (kod)^. The next dom¬ 
inant terms are the magnetic dipole, oc k^d, and electric 
quadrupole, oc kod. However, when l 2 dl = —/id/, the 
dimer only excites the modes with n being an 

even number, and the modes with n being an 

odd number. Therefore, the net electric dipole 
mode) is canceled out, and the radiation is then domi¬ 
nated by the magnetic dipole mode) and electric 

quadrupole mode) terms. 

To finalize this analysis, let us inspect the maximum 
value of the backscattering cross-section. As derived in 
Appendix A, when the response of the individual par¬ 
ticles can be approximated by their dipolar response, 
the upper bound of the backscattering cross-section of 
a nanoparticle dimer is given by 

(-k) < 1 ^ 4.5 (9) 

Therefore, the upper bound (9) represents an enhance¬ 
ment by a factor 6.25 with respect to the bound of a 
single dipolar particle: ( 1 . 5 )^/ 7 r 0.72. Thus, in theory, 

nanoparticle dimers might present significant backscat¬ 
tering enhancements as compared to individual particles. 
However, as demonstrated in [37], the upper bound of the 
backscattering cross-section for a nanoparticle system in 
which a given number of multipoles: = 2 , are excited 

is 

cTb (-k,) < ^ {Nl + 2Nj,f Ci 5.1 (10) 

Consequently, in terms of backscattering, it is concluded 
that a nanoparticle dimer cannot fully exploit the po¬ 
tential of all of the dipole and quadrupole terms simul¬ 
taneously, i.e., the aforementioned electric quadrupolar, 
magnetic dipole, TE^^^, modes cannot be 
optimally excited at the same time. This suggests that 
the scattering properties of coupled nanoparticle systems 
can be more intuitively assessed by using antenna array 
formulations, as was done with ( 2 ), instead of employ¬ 
ing the entire multipolar decomposition. Moreover, one 
could then draw upon previous superdirective antenna 
results to guide the understanding of dimer superscatter¬ 
ing behaviors. 


Z 2 I 2 = Eodl e ^ ^ Zc /i ( 12 ) 

This system of equations corresponds to the equivalent 
circuit model depicted in figure 1(b). In this circuit 
model, the source voltages are represented by the local 
electric field acting on each of the particles, integrated 
over the differential length of the current moment that 
reproduces its dipolar response, i.e., Iidl = \pi\/ joo and 
l 2 dl = \p 2 \/juj. Similarly, the current flowing in the cir¬ 
cuit represents the magnitude of these current moments, 
/i and I 2 . Therefore, Zi and Z 2 are the self-impedances 
associated with each nanoparticle current, and Zc is 
the coupling impedance, i.e., the usually written mutual 
impedances and their reciprocal relation Z 21 = Z 12 are 
replaced by it. In general, the impedance term of each 
particle is given by Zp = —r]o{kodl)‘^/ ( 67 r 6 ff^), where 
bii^ is the scattered field coefficient of the nanoparticle’s 
electric dipole (TM®^^) mode [40]. 

This simple formulation serves to elucidate the re¬ 
sponse of each individual nanoparticle, their interaction, 
and, consequently, the nanoparticle dimer. For exam¬ 
ple, in lossless particles, the resistance reduces to the 
scattering resistance of the dipole mode: Ri = R 2 = 
Edp = ( 6^)5 which is independent of the 

nanoparticle’s properties. In contrast, the reactance is 
strongly determined by the nanoparticle’s characteris¬ 
tics. For example, electrically small dielectric particles 
with radius ri and relative permittivity 5 are character¬ 
ized by the capacitive reactance X = l/(— cjC), with 
C = 47r 5o (rf/d/^) [(5 — 1) / (5 + 2)]. On the other 
hand, electrically small plasmonic nanospheres with rel¬ 
ative permittivity described by the lossless Drude model: 
6 = 1 — uOp/uo‘^^ are characterized by a series LC cir¬ 
cuit whose reactance is: X = uL + l/(— cjC), with 
L = 3 (d/^/rf) /(47r Sq cj^), ujp being the plasma fre¬ 
quency, and C = 4:7t Sq (rf /d/^ ). 

Next, the local field acting on the nanoparticle is given 
by the addition of the incident electric field plus the scat¬ 
tered electric field produced by the neighboring nanopar¬ 
ticle. The latter is described via the coupling impedance, 
Zc, which can be determined by inspecting the electric 
field excited by an EHD and dividing it by its current 
moment: 







^-3 kod 


(13) 
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For small separation distances, Zq can be approximated 
as: 


Zc = Rc + 3^c 


~^dp Rd ^ j 


1 

Zc 


(14) 


Thus, the reactance of the coupling impedance is dom¬ 
inated by a negative capacitance (positive reactance), 
with C 47r 5o (R/dl‘^. In particular, the dispersion 
profile of the coupling reactance is inverted with re¬ 
spect to that of the particle’s self-impedance. This en¬ 
sures that the structural (dimer) resonances only occur 
at discrete frequencies. Moreover, it also means that the 
resonance phenomena cannot extend over an unlimited 
bandwidth in passive particles. In addition, there is a 
small and negative resistance whose magnitude approxi¬ 
mately equals the dipolar scattering resistance. In order 
to avoid singularities it is convenient to point out that 
the absolute value of this resistance is actually slightly 
smaller. This correction value can be approximated as 
Rs ^ Rdp {kod)‘^ /5 > 0. 

The solution to the system of equations (11)-(12) pro¬ 
vides the current moments excited in the nanoparticles 


h 


Z2 + Zce-^^^^ 

Z1Z2 - zl 


. k^d 

Eodl e^ — 


(15) 


I 2 


Zl + 

^1^2 - zh 


■ kn d 

Eodl 


(16) 


It is apparent from (15), (16) that the resonance condi¬ 
tion, in which strong currents are excited in both parti¬ 
cles despite their electrically small size, is approximately 
given by X 1 X 2 = X^. In addition, the reactances re¬ 
quired to satisfy the condition C = —I 2 are given by 

V = -Xc-E (17) 


X2 = -Xc + t:X ( 18 ) 

kod 

It can be readily checked that the impedances (17), 
(18) match almost perfectly with the resonance condi¬ 
tion X 1 X 2 = X^. Therefore, it can be concluded that 
it is possible to excite very strong and opposing cur¬ 
rents in coupled nanoparticle systems. Moreover, the 
circuit model formulation provides us with a simple way 
to design the superbackscattering nanoparticle dimers: 
First, it is evident from (17), (18) that the particles must 
feature a capacitive reactance, meaning that the super¬ 
backscattering resonance is shifted toward smaller fre¬ 
quencies (longer wavelengths) with respect to the parti¬ 
cle resonance. Second, since \Xc\ ^ \Rs/kod\, it follows 
that nanoparticle dimers composed of identical particles 
actually would approximately fit to the superbakcscat- 
tering condition. Therefore, they can be considered as a 
good starting point in the design process. Finally, since 


Xi < X 2 in (17), (18), the design can be refined further 
by increasing the resonance frequency of the first parti¬ 
cle (or decreasing the resonant frequency of the second 
particle). We emphasize that this design process is quite 
general, and it can be applied to any dimer constituted 
by dipolar particles. 


IV. PLASMONIC NANOPARTICLE DIMERS 

As a first example of the use of the proposed method¬ 
ology, we study the implementation of dimers composed 
of plasmonic nanoparticles. Specifically, we assume that 
the plasmonic particles are made of aluminum-doped zinc 
oxide (AhZnO) semiconductors [53, 54], characterized by 
a relative permittivity that is described by a lossy Drude 
dispersion model: e = 1 — — jujujc)^ with plasma 

frequency uOp = 27r-2213.2 THz. For the case of simplicity, 
the lossless, ujc = 0, case is considered first, though realis¬ 
tic losses uJc = 2co’p X10“^ will be addressed later. The rel¬ 
atively low-loss plasmonic response of AhZnO is particu¬ 
larly convenient for the design of nanoparticle dimers. In 
fact, subwavelength AhZnO nanoparticles exhibit strong 
dipole resonances at near-IR and short-wavelength IR 
frequencies. Another attractive aspect is that the plasma 
frequency of AhZnO (and hence the resonance frequency 
of the nanoparticles) can be tuned by adjusting the dop¬ 
ing level [53]. 

The starting point of the design process consists of a 
dimer with identical particles of vq = 100 nm radius, sep¬ 
arated by the distance d = 3ro. Following our circuit 
model approach, figure 2 depicts the reactance of the in¬ 
dividual nanoparticles composing the dimer, Xi = X 2 , 
as well as the coupling reactance between them, Xc- As 
expected, the frequency dispersion of the reactance of 
the plasmonic nanoparticles approximates that of a se¬ 
ries LC circuit (capacitive at low frequencies, while in¬ 
ductive at frequencies above the resonance, i.e., a change 
in sign from one side of the resonance to the other). 
The resonance frequency of each individual nanopar¬ 
ticle {Xi = X 2 — 0) takes place approximately at 
A = 1537 nm, and it is denoted by the vertical black dot¬ 
ted line. On the other hand, the vertical green dashed- 
dotted line indicates the dimer resonance, Xi + Xc — 0, 
which takes place at longer wavelengths A = 1564 nm. 
The backscattering spectra of both the single nanoparti¬ 
cle and the nanoparticle dimer are depicted in figure 2(b), 
where backscattering peaks appear at the single nanopar¬ 
ticle and the nanoparticle dimer resonances, respectively. 
In addition, the backscattering cross-section peak of the 
nanoparticle dimer is 3.04, i.e., approximately 4.24 times 
larger than the backscattering bound (1) for the dipolar 
response cr 5 (k^) = (1.5)^/7r 0.72. 

This example confirms that dimers composed of iden¬ 
tical nanoparticles actually are excellent backscatterers. 
However, this 3.04 value is still below the theoretical 
bound (9). Consequently, the dimer should be optimized 
further. In this regard, the circuit model approach, i.e.. 
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FIG. 2. (a) Particle, Xi, and coupling, Xc, reactances of a 
nanoparticle dimer composed of identical lossless plasmonic 
particles of ro = 100 nm radius, separated by the distance 
d = 3ro. (b) Backscattering cross-section of the single plas¬ 
monic particle and the nanoparticle dimer. Vertical black 
dotted line indicates the single particle resonance frequency 
Vi ~ 0. Vertical green dashed dotted line indicates the dimer 
resonance: Xi H- Xc — 0. 


equations (17), (18), suggests that the backscattering re¬ 
sponse of the dimer can be enhanced by shifting the res¬ 
onance of the first particle towards higher frequencies. 
Fortunately, the resonance frequency of AhZnO nanopar¬ 
ticles can be tuned by adjusting the plasma frequency 
through the doping level. Therefore, the dimer can be 
optimized by assigning to the first particle a plasma fre¬ 
quency Upi = Up {1 Acjp), larger than the plasma fre¬ 
quency of the second particle Up. The colormap in fig¬ 
ure 3(a) represents the backscattering cross-section peak 
as a function of the inter-particle separation, d, and the 
plasma frequency deviation Aup. Being consistent with 
the circuit model, the backscattering cross-section peak 
at each separation distance is optimized by a different 
plasma frequency deviation Aup indicated by the dashed 
black line. In addition, this optimal backscattering peak 
approaches the bound (9) as the inter-particle separation 
tends to zero. 

Although our methodology identifies the optimal con¬ 
figurations that reach the theoretical bounds in the ideal 
lossless case, the dissipation losses cannot be neglected in 




A(d 

P 


(b) 

FIG. 3. Golormap of the backscattering cross-section peak of 
(a) lossless (uc = 0) and (b) lossy (uc = 2ct;pXlO“^) plasmonic 
nanoparticle dimers, as a function of inter-particle separation 
d and plasma frequency deviation Aup. The dashed black 
line indicates the position of the backscattering maximum for 
each separation distance. 


the design of plasmonic nanoparticle architectures. For 
instance, increasingly smaller separation distances be¬ 
tween the particles are inevitably associated with higher 
Q resonances, and, consequently, a higher sensitivity 
against losses. Therefore, in practice, the backscatter¬ 
ing cross-section is not optimized for the smallest sep¬ 
aration distance. In order to illustrate this fact, fig¬ 
ure 3(b) depicts the backscattering cross-section peak as 
a function of the inter-particle separation, d, and the 
plasma frequency deviation Aup, for realistic losses Uc = 
2up X 10“^. In contrast to the lossless case, the backscat¬ 
tering cross-section peak is optimized at the specific sep¬ 
aration distance d = 4.5ro, and plasma frequency devia¬ 
tion Aup = 0.025. At this optimal point, the backscat¬ 
tering cross-section approximately equals 2.8. Although 
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this value is smaller than the upper bound (9), it nev¬ 
ertheless represent a 3.9 enhancement factor with re¬ 
spect to the individual nanoparticle bound. Therefore, it 
can be concluded the superbackscattering nanoparticle 
dimers can be successfully designed with realistic plas- 
monic nanoparticles. 

In order to validate our methodology, figure 4(a) de¬ 
picts a comparison of the backscattering spectra obtained 
at the optimal point {d = 4.5ro, = 0.025) with the 

circuit model (CM) methodology, and with COMSOL 
Multiphysics 5.0 (Num.) . It is apparent from the figure 
that there is an excellent agreement between both meth¬ 
ods, validating the circuit model approach. In order to 
complete the description of the dimer superbackscatter¬ 
ing behavior, figure 4(b) depicts the electric field at the 
dimer resonance (A = 1575 nm). Note that the parti¬ 
cles appear cut in half due to the symmetry planes em¬ 
ployed in the numerical simulation. The observed elec¬ 
tric field distribution confirms the fact that the optimal 
backscattering response is based on the excitation of elec¬ 
tric dipoles of similar magnitude but opposite directions. 
Moreover, figure 4(c) represents the scattering directiv¬ 
ity pattern in the XZ- (solid line) and XY- (dashed line) 
planes at the dimer resonance (A = 1575 nm). As noted 
by the bound (1), the superbackscattering response is as¬ 
sociated with balanced forward-backward superdirective 
patterns, T)scat(~ki) T)scat(ki) ^3.4. This result illus¬ 
trates the fact that the superbackscattering response is 
constructed by the simultaneous action of the extraction 
(forward scattering) and re-radiation backward scatter¬ 
ing mechanisms. 


V. HIGH PERMITTIVITY NANOPARTICLE 
DIMERS 

Pairs of plamonic nanoparticles are but one of the 
many implementations of superbackscattering dimers. 
As a matter of fact, the dimers described in this work 
can be composed of any subwavelength particle exhibit¬ 
ing a strong dipolar resonance. As an example of this 
generality, this section addresses the design of super¬ 
backscattering dimers based on the magnetic dipole res¬ 
onances occurring in high permittivity particles. Specif¬ 
ically, we make use of particles made of Tellurium (Te). 
The choice of Te is motivated by its high refraction in¬ 
dex (e.g., riTe = ^Te ~ with ^ 5.56 at 10 /im 
[55]), and its very low losses in the thermal infrared. In 
fact, spectroscopy studies reveal that Te losses might be 
as small as 0.0004 at 10 /im [55]. However, exper¬ 

iments involving the scattering from particles fit to losses 
that are two orders of magnitude larger: n^g 0.04 at 
10 /im [56]. Thus, in order to have a complete perspective 
on the impact of losses in Te particles, low n^g 0.0004, 
high n^g 0.04, and even intermediate n^g 0.004 
losses are studied here. 

The design of the Te particle dimer can be carried 
out by using exactly the same procedure adopted for 
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FIG. 4. (a) Backscattering cross-section spectra of the opti¬ 
mized plasmonic dimer and the individual nanoparticle. Com¬ 
parison of the circuit model (CM) and COMSOL (Num.) pre¬ 
dictions. (b) Electric field distribution in the XZ plane at the 
dimer resonance (A = 1575nm). The particles appear cut in 
half, which is due to the symmetry planes employed in the nu¬ 
merical simulation, (c) Scattering directivity pattern in the 
XZ- (solid line) and XY- (dashed line) planes at the dimer 
resonance (A = 1575nm). 


the plasmonic nanoparticle dimers. In this manner, the 
design starts with a dimer composed of identical parti¬ 
cles. Even the same circuit model approach can be used, 
provided that the particles self-impedance is changed to 
Zp = —r]o{kodl)‘^/ (bTT where is the scattered 

field coefficient of the particle’s magnetic dipole 
mode [40]. In this manner, the design starts with a dimer 
composed of identical particles. The radius of each parti¬ 
cle is initially set at ri = r 2 = 0.88 /im, so that the mag- 
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netic dipole resonance of each individual particle takes 
place at A lO/iui. Next, the backscattering response 
is refined by increasing the resonance frequency of the 
first particle. In this case, this is simply accomplished by 
reducing its size, i.e., ri = 5 x r 2 , where 5 is a scaling 
factor to be optimized. 

Figures depicts the backscattering cross-section peak 
as a function of the inter-particle separation, d, and scal¬ 
ing factor s. The figure includes the results of (a) low 
0.0004, (b) intermediate 0.004, and (c) high 

n^g 0.04 losses. In the low loss case (see figure5(a)), 
the backscattering cross-section approaches the 4.5 value 
of the bound (9) as the separation distance decreases. In 
contrast, the backscattering cross-section is optimized at 
specific pairs of the separation distance and scaling fac¬ 
tor for the intermediate and high losses (see figures 5(b) 
and 5(c), respectively). In particular, the optimal pairs 
are d = 3.25 ro and s = 0.987 for intermediate losses, 
and d = 4.5 ro and s = 0.985 for high losses. In general, 
the larger the losses, the larger the optimal separation. 
Naturally, the maximal backscattering cross-section de¬ 
creases along with losses, and the intermediate and high 
loss cases feature global maxima of 3.05 and 1.2, respec¬ 
tively. These values correspond to enhancement factors 
of 4.24 and 1.67, respectively. Therefore, our analysis 
suggests that, even in the presence of moderate to high 
losses, superbackscattering dimers composed of high per¬ 
mittivity particles can be successfully implemented. This 
conclusion can be extrapolated to any frequency regime 
in which high permittivity materials are available. 

Again, the design procedure is validated with the use of 
numerical simulations. To this end, figure 6(a) depicts a 
comparison between the predictions on the backscatter¬ 
ing spectra of the optimized Te-particle dimer with in¬ 
termediate losses, obtained with the circuit model (CM) 
and COMSOL Multiphysics 5.0 (Num.). The spectra of 
the individual particle is also included for the sake of 
comparison. The figure shows that there is a good agree¬ 
ment between both approaches, validating the use of the 
CM. However, the agreement is not as excellent as it was 
for the plasmonic dimer (see figure 4(a)). The reason be¬ 
hind this slightly larger disagreement is the excitation of 
electric dipoles in the individual particles. Note that the 
CM only takes into account the magnetic dipole mode ex¬ 
cited in the Te particles, and neglects the electric dipole 
response, not necessarily small in a high permittivity par¬ 
ticle. While figure 6(a) reveals that this is actually a good 
approximation near the magnetic dipole resonance, the 
accuracy of the CM is even better for particles whose re¬ 
sponse is dominated by the electric dipole (e.g., the plas¬ 
monic particles presented in 4(a)). This outcome occurs 
because of the scattering responses omitted by the circuit 
model (i.e., the magnetic dipole and HOMs excited by 
the plasmonic particles) are significantly smaller. Conse¬ 
quently, the predictions of the CM are reasonably accu¬ 
rate in all the presented cases, enabling the intuitive and 
computationally inexpensive design of superbackscatter¬ 
ing dimers. 



Scaling 

(c) 


FIG. 5. Colormap of the backscattering cross-section peak of 
the Te particle dimers, as a function of the particle separation 
d and the scaling factor s of the first element, (a) Low ~ 
0.0004, (b) intermediate ~ 0.004, and (c) high ^ 0.04 
losses. The dashed black line indicates the position of the 
backscattering maximum for each separation distance. 

















To shed yet more light onto the polarization processes 
excited in the dimer, figure 6(b) represents the electric 
field distribution at the dimer resonance (A = 10.03 /im). 
In this case, the rotational character of the fields evi¬ 
dences the excitation of out of phase magnetic dipole mo¬ 
ments. To finalize the description, the scattering directiv¬ 
ity patterns at the XZ- (solid line) and XY- (dashed line) 
planes at the dimer resonance (A = 10.03 /im) are de¬ 
picted in figure 6(c). Again, the superbackscattering re¬ 
sponse is accompanied by superdirective scattering point¬ 
ing to both the forward (T)scat(ki) ^ 3.1) and the back¬ 
ward {Dscat{—^i) — 3.3) directions. 


VI. CONCLUSIONS 

This work has addressed the design of superbackscat¬ 
tering dimers. We analytically derived the conditions un¬ 
der which the backscattering cross-section of a nanopar¬ 
ticle dimer is maximized. Interestingly, the largest 
backscattering response is observed exactly when the 
dipole excited in the particles composing the dimer are of 
the same magnitude but opposite directions. When this 
condition is met, nanoparticle dimers can theoretically 
feature a backscattering cross-section 6.25 times larger 
than that of the individual particles. Moreover, it was 
demonstrated that nanoparticle dimers approaching this 
enhancement factor can be straightforwardly designed by 
using a simple circuit model. Naturally, the performance 
of these systems will be limited in practice by dissipa¬ 
tion damping and other practical issues. Despite this 
fact, our examples demonstrate that, even when real¬ 
istic losses are taken into account, both plasmonic and 
high-permittivity particle dimers can exhibit fourfold en¬ 
hancements of the backscattering cross-section with re¬ 
spect to the individual particles. These results encourage 
us to believe that the development of superbackscatter¬ 
ing nanparticles dimers is possible. Due to the comple¬ 
mentary interest of backscattering measurements, we are 
confident that these nanoparticle architectures will be of 
great interest for variety of spectroscopy, communication, 
remote sensing, manipulation and/or imaging systems. 


Appendix A: Derivation of the optimal excitation 
coefficients 



A (iLinn) 
(a) 



(b) 

90 4 



(c) 

FIG. 6. (a) Backscattering cross-section spectrum of the op¬ 
timal (d = 3.25ro, s = 0.987) Te-particle dimer with interme¬ 
diate losses {riTe = 5.56 — j0.004). Comparison between the 
predictions of the circuit model (CM) and COMSOL (Num.). 
(b) Electric field distribution in the XZ plane at the dimer 
resonance (A = 10.03/xm). The particles appear cut in half 
due to the symmetry planes employed in the numerical sim¬ 
ulation. (c) Scattering directivity pattern in the XZ- (solid 
line) and XY- (dashed line) planes at the dimer resonace. 


This appendix summarizes the maximization process 
that identifies the optimal excitation coefficients and up¬ 
per bounds of the backscattering and forward scattering 
cross-sections of electrically small nanoparticle dimers. 
To this end, and in contrast with figure 1(a), here we as¬ 
sume that the particles composing the dimer are located 
along the x-axis, and the incident field is z-linearly po¬ 
larized. This new coordinate system is particularly con- 
vement because the scattered electric field in the far zone 
is 6 polarized. Specifically, the scattering pattern of the 


dimer in this new coordinate system is given by 

f (f) = 0 sin6» (Al) 

Note also that the scattering directivity (3) can be writ¬ 
ten alternatively as: [44, 46] 


-^scat 


3 0 

-sin 0 -^^- 

2 \Aif+ \A2f+3HRe{AiA*2} 


(A2) 
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with H = H (kod) defined as 


H{kod) 


sin {kod) 

knd 


{kody 


+ ,a3) 

{kodf 


Next, in order to study electrically small dimers we 
observe the limit 


lim H (kod) 

/cq d —^0 



(A4) 


In this manner, the scattering directivity of an electri¬ 
cally small {kod 0) dimer approximately reduces to 


n |Ai + ^2 (1+jA:o(ff-x)! 

2 \A,+A2f-l{kodfRe{A,A*,} 

(A5) 

If we take the kod ^ 0 limit in (A5) with constant coef¬ 
ficients, then the scattering directivity reduces to a dipo¬ 
lar pattern I)scat(?) — (3/2)sin^6>. Therefore, in order 
to get a response different from a dipole pattern from an 
electrically small dimer, the excitation coefficients must 
increase as the separation distance decreases. The coef¬ 
ficients that satisfy this condition, while keeping a finite 
scattering pattern, can in general be written as: 


A-i 


^0 A - ( ^ ^ 

— w—^2 — ao I T + ai 

kod \kod 


(A6) 


By using these coefficients, the scattering directivity of 
an electrically small dimer is given by 


-^scat {^) 



\ai +jr-xl 




(A7) 


This equation can be maximized at any specific direction 
by using standard techniques. For example, it is clear 
from (A7) that the forward scattering directivity, r • x = 
1, is maximized with an imaginary coefficient ai = jci. 
Maximization with respect to ci leads to the solution 


Cl =2/5. By using this value, it is found that the forward 
scattering directivity of an electrically small nanoparticle 
dimer is upper bounded by 

-Dscat (x) < ^ = 5.25 (A8) 

This analytical result is consistent with the numerical 
optimizations carried out in [45, 46]. 

On the other hand, the forward-backward directivity 
product can be written as: 


^ ^ + — 4Im{ai}^ 

^scat (x) T^scat (~x) = T ~ r2 

' + i) 

(A9) 

It is apparent from (A9) that the forward-backward di¬ 
rectivity product is maximized for ai being a real num¬ 
ber. The maximization for such a real number results 
in the optimal value ai = 0, i.e., the forward-backward 
directivity product is maximized when the excitation 
coefficients satisfy A 2 = —Ai = ao/kod. This condi¬ 
tion physically corresponds to the excitation of dipoles 
with the same magnitude but opposite directions. At 
this maximal point we find the upper bound of the 
forward-backward directivity product of an electrically 
small nanoparticle dimer to be: 


-^scat (x) -^scat ( x) ^ 



14.06 


(AlO) 


Consequently, the upper bound of the backscattering 
cross-section of an electrically small nanoparticle dimer 
can be found by introducing (AlO) into (1). It is: 

CTb (-ki) < - (^^p^ = 4.48 (All) 
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